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Abstract 



This paper presents a new algorithm for the convex hull problem, which is based on a 
reduction to a combinatorial decision problem CompletenessC, which in turn can 
be solved by a simplicial homology computation. Like other convex hull algorithms, 
our algorithm is polynomial (in the size of input plus output) for simplicial or simple 
j> ■ input. We show that the "no"-case of CompletenessC has a certificate that can 

O : be checked in polynomial time (if integrity of the input is guaranteed). 

! 1 Introduction 

CO 

Every convex polytope P C M. d can be described as the convex hull of a finite set V of 
points or as the (bounded) set of solutions of a finite system Tt of linear equations and 
inequalities [23, Lect. 1]. In view of the fundamental role that polytopes play in Euclidean 
geometry and hence for any type of geometric computing, the conversion between the two 
types of representations, known as the convex hull problem, is of key interest. It splits 
^ ■ into two separate tasks. 

The first task is the facet enumeration problem: Given a finite set of points V C M. d , 
determine the combinatorial structure of its boundary. For this one does not want to 
explicitly enumerate all the faces (the intersections of P with supporting hyperp lanes) , 
but one wants sparser date, namely to compute a minimal representation of the convex 
hull conv("P) in terms of equations and (facet-defining) inequalities. Here the equations 
should describe the affme hull aff(P), while the additional inequalities correspond to the 
facets (faces of codimension 1) of P. If P is full-dimensional in M. d , then the facet-defining 
inequalities are unique up to scaling. 

The second task is the vertex enumeration problem: Given a finite system 7i of linear 
(equations and) inequalities, and provided that the set of solutions P = f] H is bounded, 
compute the minimal set of points V whose convex hull is P. This minimal set is unique; 
it consists of the vertices (O-dimensional faces) of P. 

*An extended abstract version of this paper, "Polytope verification by homology verification," has 
appeared in the Proceedings of EuroCG, Berlin, March 26-28, 2001, pp. 142-145. 
+ Partially supported by Deutsche Forschungs-Gemeinschaft (DFG) 
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The two tasks are dual to each other, via cone polarity. Thus if an LP-type oracle (an 
algorithm which for a system of inequalities computes a solution, or for a set of points 
computes a separating hyperplane, cf. [14]) is available, every algorithm for the facet 
enumeration problem can also be used for vertex enumeration, and vice versa. 

Despite the great interest in the convex hull problem, and despite the fact that a 
number of different strategies and algorithms have been explored, implemented and an- 
alyzed in detail (see [10], as well as Avis [1] [2], Fukuda [11] and Gawrilow & Joswig 
[12] [13] for implementations), the problem can be considered "solved" neither in the- 
ory, nor in practice. If the dimension d is fixed, Chazelle's celebrated algorithm [8] gives 
an asymptotically worst-case optimal (polynomial time) theoretical solution. Its opti- 
mality is based on McMullen's "Upper Bound Theorem" [20] on the maximal number 
of facets for a d-polytope with n vertices. However, for any given convex hull problem, 
the output may be small, but it may also be much larger than the input — indeed, it 
may be of exponential size, if the dimension is not fixed. This is very relevant, since 
high-dimensional computations occur in a variety of important applications. Thus one 
is asking for a convex hull algorithm whose running time is bounded by a polynomial in 
the size of "input plus output"? Such an algorithm would be called output- sensitive. The 
analysis by Avis, Bremner and Seidel [3] shows that, unfortunately, none of the known 
types of convex hull algorithms is output-sensitive. These can roughly be categorized as 
follows: Incremental and triangulation producing (e.g., Chazelle's method), incremental 
without triangulations (e.g., Fourier-Motzkin elimination [23, Lect. 1]), non-incremental 
(e.g., reverse search [4]). Note that, by a result of Bremner [6], only non-incremental 
methods can possibly be output-sensitive. 

The purpose of this paper is to describe a new (non-incremental) convex hull algo- 
rithm, based on a completely different principle. To this end, we first present a (folklore) 
polynomial reduction of FacetEnumeration to the decision problem PolytopeVeri- 
fication. Then we further reduce to the Completeness problem: Is a given description 
of a rf-polytope by some of its vertices and some of its facets complete, that is, are we 
given all the vertices and all the facets? Looking at the convex hull problem via its 
reduction to Polytope Verification or Completeness automatically reveals its in- 
herent self-dual structure. It is an interesting feature that the Completeness problem 
can be posed both with geometric input data and as an entirely combinatorial problem 
CompletenessC, where only the incidences between vertices and facets are given. 

Let us just mention here one recent occurrence of the combinatorial completeness prob- 
lem: McCarthy et al. [19] describe a situation where one wants to know whether a given 
inequality description for a polytope is complete. Moreover, the vertex coordinates in 
some of their problems are necessarily non-rational, so any coordinate-free/combinatorial 
approach is welcome. Unfortunately, the most interesting case left "open" by McCarthy 
et al. (the convex hull of the matrices corresponding to the Coxeter group H±) is a poly- 
tope completeness problem in dimension d = 16 with 14,400 vertices: From this data our 
method generates gigantic boundary matrices that are plainly too large to process. 

Also we have been informed by Samuel Fiorini (email, January 2002) that he has 
successfully used a certificate for the "no"-case of CompletenessC that is similar to the 
one that we describe in Section 6. 

Our main contribution is an algorithm to attack the combinatorial CompletenessC 
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problem via deciding whether a certain simplicial homology group of a certain abstract 
simplicial complex vanishes or not. Moreover, we present a polynomially checkable certifi- 
cate for non-completeness, provided that the input is valid. For the geometric version the 
validity of the input can be checked easily. Unfortunately, the complexity status for the 
homology computation problem is open. The best currently available strategy to decide 
non-triviality of a (rational) homology group in question seems to be to compute boundary 
matrices and perform Gauss elimination. Since the boundary matrices in our algorithm 
can be exponentially large, we do not obtain an output-sensitive method. However, like 
other methods (e.g., Avis' and Fukuda's reverse search [4] or Seidel's gift-wrapping algo- 
rithm [22]) our algorithm is output-sensitive in the case of simplicial polytopes. 

2 FacetEnumeration via PolytopeVerification 

We start with a more formal description of the facet enumeration problem: 
FacetEnumeration (e, V): 

Input : integer e > 0; finite set of points V C 1R 6 . 

Output: minimal description of conv(P) in terms of equations (for the afline hull of V) 
and inequalities (one for each facet of conv(P)) 

It is known, cf. Avis, Bremner & Seidel [3], Fukuda [10, Node 21], and Kaibel & Pfetsch [17, 
Problems 1-3], that FacetEnumeration has a polynomial reduction to the polytope 
verification problem: 

POLYTOPEVERIFICATION(e, V, Ti) : 

Input : integer e > 0; finite set of points PcK e ; finite set H of closed halfspaces in W 
Output: answer yes/no to the question whether conv(V) — f\H 

Freund and Orlin could show that a related problem, to decide whether f]7i C 
conv('P), is co-NP-complete [9]. 

3 PolytopeVerification via CompletenessG 

Assuming that an LP-type oracle is available, the PolytopeVerification problem is 
polynomially equivalent to the following geometric polytope completeness problem: 
CompletenessG (d, V, F): 

Input : integer d > 0; finite set of points V C M d ; finite set F of closed halfspaces in R d , 
such that 

• P := conv(V) is contained in Q := f]^ 

• dim P = dim Q = d 

• every dgV defines a vertex of Q 

• every F e T defines a facet of P 
Output: answer yes/no to the question whether P = Q 

As in the case of PolytopeVerification, the roles of vertices and facets are inter- 
changeable for CompletenessG. 

We sketch the reduction of PolytopeVerification to CompletenessG. Given 
any input (e,V,H) for POLYTOPEVERIFICATION, set P := conv(V) and Q := f]H. 
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Employ Gaussian elimination to determine dim P. Verify whether all the inequalities 
in H are valid for P; if this is not the case, then P % Q, so we output no; otherwise 
P C Q is established. Now extract the set 7i' of all halfspaces from 7i for which P lies 
in the bounding hyperplane, that is, all those inequalities which are tight on aff P. An 
LP-type oracle is sufficient, but also needed [14], to check whether [^H' — aff P; if this 
is not the case, then we know that dimQ > dimP, so we can output no. Otherwise we 
proceed by restricting the input to aff P, that is, we deal with the situation where P is 
full-dimensional. 

Now remove from 7i all the halfspaces which do not determine facets of P; this may 
be done using Gaussian elimination. (In the case P = Q, this removal does not change Q; 
in the case P C Q, it may enlarge Q.) Similarly, we now remove from V all those points 
which do not arise as intersections of some bounding hyperplanes of halfspaces in H; again 
this may be done via Gaussian elimination. (In the case of P = Q, this removal does not 
change P; in the case P C Q, we may loose vertices of P, thus making P smaller.) 

Now we have prepared our input for CompletenessG. Indeed, the first two conditions 
on the input are satisfied, the other two are easily checked: If one of them fails, then output 
the answer no. □ 



4 CompletenessG via CompletenessC 

The incidence matrix of a polytope P with vertex set V and facet set T is defined to be 
the matrix 

Ip := (iFv)F<=J r ,v£V 

where ip v — 1 if vertex v lies on the facet F (that is, if v e F), and ip v = means that 
v F. This matrix is well-defined up to permutation of rows and of columns, which 
corresponds to reordering V and T . A minor of a matrix will refer to any submatrix 
obtained by possibly removing rows and/or columns. A minor J of the incidence matrix Ip 
is complete if J = Ip. Thus we arrive at the combinatorial polytope completeness problem: 
COMPLETENESSC(<i, J): 

Input : integer d > 0; incidence matrix minor J of a d-polytope 
Output: answer yes/no to the question whether J is complete 

It is not obvious that this problem is well defined. However, from Theorem 5.1 below it 
follows that there are no two <i-polytopes P and P' such that a 0/1-matrix J is both a 
complete incidence matrix for P and an incomplete minor of an incidence matrix for P' . 
(See also the related discussion in [16].) It is clear that CompletenessG has a polynomial 
reduction to CompletenessC. 

It is essential to have the dimension among the input parameters of CompletenessC. 
This is demonstrated by the following example [23, p. 71]: 

1 1 \ 
110 
110 1 
10 110 
11110 
1111/ 
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We can identify V = {1, 2, . . . , 8} and T = {1234, 1278, 1458, 2367, 3456, 5678} with the 
sets of vertices and facets, respectively, of a 3- dimensional cube (in a suitable "Klee-Minty" 
vertex numbering; see Figure 1(b) below). Consequently, COMPLETENESSC(3, Jkm) — 
yes. But we can also identify V with the vertices of a cyclic 4-polytope 6*4(8). Then 
each element in T corresponds to a facet of 6*4(8), according to Gale's evenness criterion. 
Hence CompletenessC(4, J km ) = no, since C 4 (8) has 20 facets. 

A more generic class of examples for which the dimension information is needed arises 
from the prism construction: Let P be an arbitrary <i-polytope and P' = P x [0, 1] the 
prism over P. The facets of P' are P x {0}, P x {1}, and the products of facets of P with 
the interval [0, 1]. Call the latter facets of P' vertical, and let Jp be an incidence matrix 
of P. We have CompletenessC(o?, J P ) = yes. On the other hand J P is also a minor of 
an incidence matrix of P', which corresponds to the vertical facets and, say, the vertices 
in the bottom facet P x {0}. Therefore, CoMPLETENESsC(<i + 1, Jp) = no. 

5 CompletenessC via simplicial homology 

We will point out that CompletenessC has a topological core. The reader is referred 
to Bjorner [5] for a survey of topological combinatorics tools, and to Munkres [21] for 
a presentation of simplicial homology. In the following we will use reduced simplicial 
homology with coefficients in Z 2 . One could use any other commutative coefficient ring 
with unit, but Z 2 is the natural choice in terms of efficiency and simplicity. We choose 
non-reduced homology to simplify notation for the trivial case d = 1. 

Let J G {0, l} jrxV be an incidence matrix minor of some polytope P with vertex 
set V D V and facet set T' D T . Thus the columns of J are in bijection with a (partial) 
vertex set V of P. Each row of J is the characteristic vector of a subset of rows, i.e., of a 
subset of V. Thus in the following we interpret J as a combinatorial encoding of a system 
T of (not necessarily distinct) subsets of V, and with slight abuse of notation we write 
T C 2 V . The crosscut complex of J is the simplicial complex 

T(J) := (V, \J{2 F :FEF}), 

the simplicial complex of all sets of vertices that are contained in some facet in T . 

Theorem 5.1. The incidence matrix minor J G {0, l} jrxV of a d-polytope is complete if 
and only if Hd-i(T( J); Z 2 ) 7^ 0. 

Proof. The set 

n(P, J) : = (J conv{v G V : v G F} C dP 

is a compact subset of the boundary of P: For every "given" facet F of P, it contains 
the convex hull of all "given" vertices. Thus H(P, J) is a polyhedral complex, called a 
partial polytope, covered by its convex (and hence contractible) cells conv{t> G V : v G F}. 
According to the nerve theorem [5], the crosscut complex T(J) has the same homotopy 
type as the set IT(P, J). In particular, the homology of the set T1(P, J) and of the crosscut 
complex coincide. For an example of the crosscut complex of a partial polytope see 
Figure 1(a). 
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In the yes case, if the sets of vertices and facets both are complete, n(P, J) is the 
complete boundary of P, homeomorphic to S^ -1 , so we have P rf _!(r(jF); Z 2 ) = Z 2 . 

In the no case, if the vertex or the facet list is incomplete, then II(P, J) is a proper 
subset of dP, which is a subcomplex of a suitable triangulation of dP, so it cannot have 
(d — l)-dimensional homology. □ 

The complexity status of the problem to compute the rank of an arbitrary homology 
group, or even to decide whether a certain homology group vanishes, seems to be open; see 
Kaibel & Pfetsch [17, Problem 33]. Thus currently our best option is based on explicitly 
computing simplicial homology via boundary matrices, as in Algorithm A. 



Algorithm A: COMPLETENESS VlAHOMOLOGY((i, J) 

Input : integer d > 0; an incidence matrix minor J of a <i-polytope 
Output: answer yes/no to the question whether J is complete 

Generate Z 2 -boundary matrices and da-i for T(J) 
if dim Z2 ker c^-i > rank^ 2 <9 d then 

j return yes 
else 

L return no 



To estimate the costs of this computation, suppose that n — |V|, m — and that 
the maximum cardinality of any facet equals s. Thus J G {0, i} mx ™ ) and every row of J 
contains at most s ones. Then the size of the relevant boundary matrices is bounded from 
above by x ( s d ) m an d Q) m x id-i) 171 ' respectively. We use Gaussian elimination 

over Z 2 to compute the rank and the corank, respectively. 

Corollary 5.2. The algorithm COMPLETENESS VlAHOMOLOGY(<i, J) has a polynomial 
running time if s is bounded by d + c, for an absolute constant c > 0. 

The latter case is, in fact, interesting: A d-polytope is simplicial if each proper face is a 
simplex or, equivalently, each facet contains exactly d vertices. We infer that the running 
time of Completeness ViaHomology for simplicial polytopes is bounded by 0(dm 3 ). 

It has been observed by Bremner, Fukuda & Marzetta [7] that FacetEnumeration 
for a polytope P is polynomially equivalent to FacetEnumeration for the dual poly- 
tope P*. Using our techniques, a similar result can be obtained directly. If I is an 
incidence matrix for P, then the transposed matrix J tr is an incidence matrix for P*. Any 
minor J of / is complete if and only if its transpose is a complete minor of P r . This leads 
to the following modification of our algorithm. While s was defined above as the maximal 
row size of the input incidence matrix minor, define 

s' := min{ maximal row size, maximal column size}. 

Thus we modify our algorithm: It should first compare the sizes of the primal and the dual 
problem, and then perform the (reduced) homology computation for the smaller problem. 
The modified algorithm COMPLETENESSVlAHOMOLOGY(d, J) has polynomial running 
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time if s' is bounded by u d plus a constant." In particular, this yields an 0(d(n + m) 3 )- 
algorithm for the CompletenessC problem specialized to polytopes which are simplicial 
or simple, that is, dual to a simplicial polytope. 

We note, however, that these running times are neither optimal nor the best available: 
The reverse search algorithm of Avis and Fukuda [4] computes the convex hull (and 
thereby solves CompletenessG) of a simplicial polytope in 0(dnm) steps. 



6 A Certificate for Incompleteness 

Let P be a rf-polytope with ordered vertex set V = {v\, . . . ,v n } and facet set T' . In- 
ductively, define a sequence A , . . . , A ra of polytopal subdivisions of the boundary com- 
plex dP: Set A := dP. In order to obtain A^ replace each facet F of A i i c _ 1 which con- 
tains Vk by the set of cones with apex over those facets of F which do not contain v^. 
The final subdivision is a triangulation A(P) := A m of dP, the pulling triangulation [18] 
with respect to the chosen ordering of V. For an example of a pulling triangulation see 
Figure 1(b). 





(a) The (3-dimensional) crosscut complex of 
some partial 3-cube C. The two quadrangle 
faces of C yield tetrahedra in T(C), which 
are displayed almost flat. 



(b) The pulling triangulation of the bound- 
ary of a 3-cube with respect to a "Klee- 
Minty" vertex ordering. The facet {1,7,8} 
of the triangulation corresponds to the flag 
{8} c {7, 8} c {1, 2, 7, 8} of the cube. 



Figure 1: Crosscut complex and pulling triangulation. 



The pulling triangulation of dP has several nice properties (not shared, for example, by 
the "placing triangulation") which may be exploited for our purposes. First, its combina- 
torics is determined by the combinatorics of P; see below. Furthermore, if we use a linear 
ordering of the vertex set V in which the vertices in V come first, then the corresponding 
pulling triangulation of the boundary of P contains II(P, J) as a subcomplex. 

Let us now identify the vertex set V' with the set [n] = {1, . . . , n} and each facet F e 
T' with the subset of [n] that corresponds to the vertices contained in F. Thus any 
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triangulation of dP is encoded by a collection of ci-subsets of [n], that is, to a subset of 
. We write {v i, . . . , v d } < for a <i-subset of [n] with vi < v 2 < • • • < v d . 

Lemma 6.1. Let P be a d-polytope whose vertex set is labeled by [n\. 

Then a set {vi, . . . ,Vd}< G (^') corresponds to a facet of the pulling triangulation of dP 

(with respect to the chosen vertex labeling) if and only if there is a complete flag of faces 

c Go c d c ... c G d -i c P, 

such that Vi is the smallest vertex in Gd-i for 1 < i < d, that is, if there are facets 
Fi, . . . ,Fd of P such that 

Vi = min(Fi n . . . D Fi) 

for 1 < i < d. 

Proof. Every pulling facet {i>i, . . . , Vd}< lies in a facet F 1 = G d -\ of P, with v\ = min G d -i- 
It is a cone with apex V\ and base G d -2 C G d -i- The existence of the rest of the maximal 
flag (Gi) o<i<d follows recursively. Given the flag, the existence of the facets F 1 , . . . , F d 
follows [23, Lect. 2]. Given a complete flag, the corresponding sequence of facets Fi is 
uniquely determined if P is simple, but not in general. □ 

If we have an arbitrary incidence matrix minor J of a d-polytope P, then we can read 
the combinatorial characterization of the pulling triangulation from Lemma 6.1 as the 
definition of a complex that coincides with the pulling triangulation of dP in case J is 
complete, but is well-defined in general: 

Definition 6.2. Given an integer d > and a 0/1-matrix J G {0,l} mxn , which we 
interpret as the incidence matrix of a set system T C 2^, the pulling complex old and J is 

A(d, J) := • • • , v d }<: G (^') : there are F 1? . . . , F d G T such that 

Vi = min(Fi n . . . n Fi) for 1 < i < dj. 

Lemma 6.3. Let P be a d-dimensional polytope with vertex set V' and facet set T' , and 
let J be a incidence matrix minor corresponding to subsets V C V' and T C T' . Let 
P C P be the convex hull of the vertices in V. Fix a linear ordering on the vertex set V' 
such that the vertices in V come first. 

Then the simplicial complex A(d, J) is a subcomplex of A(P) as well as of A(P). 
In particular, A(d, J) is a proper subcomplex of A(P), unless the minor J is complete, 
J = Ip. In the incomplete case A(d, J) may even be empty. 

Proof. Let {v i, . . . , v d }< G A(d, J), then there are Fi, . . . , F d G T such that Vi = min(Fi D 
. .. fl Fi). Now since J is an incidence matrix minor of P, there are facets Fi D F 
of P, and by the assumption on the vertex ordering the vertices in F come first, so 
min(Fi H . . . n Fi) = min(Fi D . . . D Fi), which yields {f i, . . . , v d } < G A(P). 

Now P = conv(V), and the Fi = Fi fl V are vertex sets of faces (not necessarily facets) 
of P. If the vertices Vi = min(Fi fl . . . fl F) are distinct, then the faces F± fl . . . fl F form a 
complete flag in the face lattice of P, and thus {v±, . . . , v d } < G A(P), by Lemma 6.1. □ 
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In particular, A(d, J) triangulates a subset of the complex IT(P, J) that appears in the 
proof of Theorem 5.1. 

Now we present a polynomially-checkable certificate for the case that J is incomplete. 
Note, however, that this result does not prove that CompletenessC is in co-NP: We 
are not able to check (in polynomial time) whether the input is valid, that is, whether J 
is actually an incidence matrix minor of some <i-polytope. 

Theorem 6.4. Any no instance of the problem CompletenessC (rf, J) has a certificate 
that can be verified in polynomial time. 

Proof. The minor J is incomplete if and only if the pulling complex A(d, J) is not a 
complete triangulation of a rf-polytope boundary. Two cases arise. The first one is if 
A(d, J) = 0, in which case Algorithm B described below will certify in polynomial time 
that J is not complete. 

The second case is if A(d, J) is non-empty but incomplete. In this case (since the 
dual graph of the pulling triangulation A(P) is connected) there is a facet {v i, . . . , v^} G 
A(d, J) together with an index % such that there is no second facet of A(d, J) that contains 
{v i, . . . , v d } \ {vi}. In this situation our certificate is the set . . . , v d } \ {vi}. Calling 
IsPullingFacet for every rf-subset of [n] which contains the certificate, this certificate 
can be verified in polynomial time, since there are n — d + 1 of these subsets. □ 

Now we proceed by describing the two subroutines needed for Theorem 6.4. The first 
one is Algorithm B: Given an incidence matrix minor J it either finds a facet of A(d, J) 
in polynomial time or it detects that J is incomplete. The correctness follows from 
Lemma 6.1. Our specific formulation of the algorithm produces a pulling triangulation 
facet which does not contain 1: This restriction does not hurt, since A(d, J) must contain 
such a facet if J is complete. 



Algorithm B: FindPullingFacet(c?, J) 

Input : incidence matrix minor J G {0, l} mxn of a rf-polytope; d-tuple 

K...,M<g( [ J) 

Output: a facet {v i, . . . , Vd} G A(d, J), or incomplete 

S<-[n] 

for % <— 1 to d do 

Fi <— any FeJ such that min S £ F, F n S ^ 0, and |F n S\ is maximal 
if no such facet exists then 
|_ return incomplete 

S <- S n F t 
_ Vi min S 

return {v u . . . ,v d } < 



Our second subroutine, Algorithm C, checks whether a given set of d vertices is a facet 
of the pulling complex A(d, J) or not. Its correctness again follows from the characteri- 
zation in Lemma 6.1. Its running time is bounded by 0(d(n + to)). 
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Algorithm C: IsPullingFacet^, J, {v ± , . . . , Ud}<) 



Input : (d, J) as above 

Output: answer yes/no to the question whether {vi, . . . , v d } G A(d, J) 
for i <— d downto 1 do 

|_ compute the set Ti of all facets (i.e., rows of J) that contain {v iy . . . , i^} 

for i <— 1 to d do 

Fj <— any F G ^ with ^ = min(Fi Pi ... Pi Fj_i fl F) 
if no swc/i F exists then 
L return no 

return yes 



We close our discussion with a pointer to a specific special case: It would be interesting 
to know whether Completeness (<i, J) has a polynomial time solution for the very special 
case where J has all columns and lacks at most one row. 

Acknowledgements 

We are grateful to Volker Kaibel, Marc E. Pfetsch and Mark de Longueville for helpful 
comments. Moreover, the first author is indebted to Giinter Rote for an enlightening 
discussion on the subject. 



References 

D. Avis: Irs — a C implementation of the reverse search vertex enumeration algorithm, 
version 4.1. http : //cgm. cs .mcgill . ca/~avis/C/lrs .html 

: Irs: A revised implementation of the reverse search vertex enumeration algorithm, 

in: "Polytopes — Combinatorics and Computation" (G. Kalai & G. M. Ziegler, eds.), 
DMV-Seminars Vol. 29, Birkhauser Basel 2000, pp. 177-198. 

D. Avis, D. Bremner & R. Seidel: How good are convex hull algorithms?, Comput. Geom. 
7 (1997), 265-301. 

D. Avis & K. Fukuda: A pivoting algorithm for convex hulls and vertex enumeration of 
arrangements and polyhedra, Discrete Comput. Geom. 8 (1992), 295-313. 

A. Bjorner: Topological methods, in: "Handbook of Combinatorics" (R. L. Graham, M. 
Grotschel & L. Lovasz, eds.), Vol. II, North-Holland, Amsterdam 1995, pp. 1819-1872. 
D. Bremner: Incremental convex hull algorithms are not output sensitive, Discrete Comput. 
Geom. 21 (199), 57-68. 

D. Bremner, K. Fukuda & A. Marzetta: Primal-dual methods for vertex and facet enumer- 
ation, Discrete Comput. Geom. 20 (1998), 333-357. 

B. Chazelle: An optimal convex hull algorithm in any fixed dimension, Discrete Comput. 
Geom. 10 (1993), 377-409. 

R. M. Freund & J. B. Orlin: On the complexity of four polyhedral set containment problems, 
Math. Programming 33 (1985), 139-145. 



10 



K. Fukuda: Frequently asked questions in polyhedral computation, Oct. 2000. 
http : //www. cs .mcgill . ca/~f ukuda/ sof t/polyf aq/polyf aq.html 

: CDD — A C-implementation of the double description method; cddlib, version 

0.92b. http : //www. cs .mcgill . ca/~fukuda/sof t/cdd_home/cdd.html 
E. Gawrilow & M. Joswig: polymake, version 1.5.1: a software package for analyzing convex 
polytopes, 1997-2003. http://www.math.tu-berlin.de/polymake 

: polymake: an approach to modular software design in computational geometry, in: 

Proc. 17th Annual Symp. Computational Geometry, ACM, 2001, June 3-5, 2001, Medford, 
MA, pp. 222-231. 

M. Grotschel, L. Lovasz & A. Schrijver: Geometric Algorithms and Combinatorial Opti- 
mization, Algorithms and Combinatorics, Vol. 2, Springer- Verlag, Berlin Heidelberg, 1988; 
second edition 1994. 

M. Joswig: Beneath- and-beyond revisited, to appear in "Algebra, Geometry, and Software 
Systems" (M. Joswig &: N. Takayama, eds.), Springer- Verlag, Heidelberg 2003, pp. 1-21. 
M. Joswig, V. Kaibel, M. E. Pfetsch & G. M. Ziegler: Vertex-facet incidences of unbounded 
polyhedra, Adv. Geometry 1 (2001), 23-36. 

V. Kaibel & M. E. Pfetsch: Some algorithmic problems in polytope theory, in: "Algebra, 
Geometry, and Software Systems" (M. Joswig & N. Takayama, eds.), Springer- Verlag, Hei- 
delberg 2003, pp. 23-48. 

C. W. Lee: Subdivisions and triangulations of polytopes, in: "Handbook of Discrete and 
Computational Geometry" (J. E. Goodman h J. O'Rourke, eds.), CRC Press, Boca Raton 
1997, pp. 271-290. 

N. McCarthy, D. Ogilvie, I. Spitkovsky & N. Zobin: Birkhoff's theorem and convex hulls of 
Coxeter groups, Linear Algebra Appl., 347 (2002), 219-231. 

P. McMullen: The maximum numbers of faces of a convex polytope, Mathematika 17 (1970), 
179-184. 

J. R. Munkres: Elements of Algebraic Topology, Addison- Wesley, Menlo Park, CA 1984. 
R. Seidel: Small- dimensional linear programming and convex hulls made easy, Discrete 
Comput. Geom. 6 (1991), 423-434. 
[23] G. M. Ziegler: Lectures on Polytopes, Graduate Texts in Math., Vol. 152, Springer- Verlag, 
New York 1995; revised 1998. 



11 



